The urban economy as a scale-free network 
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We present empirical evidence that land values are scale-free and introduce a network model that 
reproduces the observations. The network approach to urban modelling is based on the assumption 
that the market dynamics that generates land values can be represented as a growing scale-free 
network. Our results suggest that the network properties of trade between specialized activities 
causes land values, and likely also other observables such as population, to be power law distributed. 
In addition to being an attractive avenue for further analytical inquiry, the network representation 
is also applicable to empirical data and is thereby attractive for predictive modelling. 
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I. INTRODUCTION 

The Zipf rank-size law for city sizes is one of the most 
widely known power laws in science 0. It is also but one 
out of many similar power laws from systems in biology, 
economy and society. We continue this research by pre- 
senting empirical evidence that land values are scale-free. 
The data we use is based on a database delivered by Swe- 
den Statistics that covers estimations of the market value 
of all land in Sweden (2.9 million data points). 

Although power laws are common they are not eas- 
ily reconstructed from realistic underlaying dynamics. 
Their ubiquity suggests that they could be caused by 
some general systemic property common to a range of 
systems. Recent research suggests that many empirically 
observed power laws may be due to fundamental prop- 
erties of these systems viewed as networks of interacting 

nodesHHH. 

We investigate the mechanisms causing 
land values to follow these statistics and present a net- 
work model that reproduces the empirical results. The 
model is based on basic definitions of city formation in 
urban economics theory |5j- 

According to urban economics theory, the formation 
of modern cities is primarily caused by the advantages 
of trade between specialized producers. The exchange of 
goods and services between localized and largely immo- 
bile activities in trade economies makes a network rep- 
resentation natural: the nodes are units of land and the 



edges represent the exchange of goods and services be- 
tween them. 

It has been shown that the node degrees of a certain 
class of growing networks are power law distributed0,|(|. 
This class of networks is important because their growth 
mechanisms can be mapped to the microscopic dynam- 
ics of several real- world systems. We demonstrate how 
trade in an urban system can be represented as a scale 
free network and that, as a consequence, land values can 
be expected to follow the same distribution. We also ver- 
ify that the model retains these properties when spatial 
constraints are taken into account. We do this by us- 
ing a spatial network model to reproduce the empirically 
observed distribution of land values. 

The network approach solves a fundamental issue in 
the problem of modelling urban systems by representing 
the system at the level of the underlaying market struc- 
ture. This allows us to produce prices in units of currency 
rather than undefined and subjective fitness measures. 
It thereby opens up doors for several extensions of the 
scope of the model and provides a natural interface for 
integrating it with other models. 

When we refer to an urban system we do not necessar- 
ily refer to individual cities but rather to systems of spe- 
cialized trading activities. To clarify further, our use of 
the term activity refers to trade gains in units of currency 
per unit area and unit time. Activities can be resolved 
to any resolution down to individual transactions. 
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FIG. 1: The nodes are non-overlapping areas of land. Areas 
with no activity are not shown in the figure. 

II. THE URBAN ECONOMY AS A NETWORK 

In section Til Al we define a non-spatial model of urban 
economic growth and in section III Bl we extend it to a 
spatial model where growth is mapped to a 2D surface. 
In section III CI we motivate the model ontology and the 
basic assumptions on which we have based the model. In 
particular we discuss the connection between node de- 
gree and land value and how the urban system meets the 
criteria for being a scale-free network. 



where n t is the number of developed nodes at 
time t and = 1 if node i is developed and 

5^ =0 otherwise. 

The second end-point is selected preferentially 
among developed nodes. Preferential selection cor- 
responds to the uniform selection of an activity in 
the system and the subsequent location of its node. 
It was defined by Barabasi et al as 



1 Ej xj ' 

where is the probability of node i to be attached 
to a new edge and Xi is the degree of node i. 

2. With probability q 2 we add edges between m pairs 
of nodes that are both selected preferentially ac- 
cording to Eq. J2J). 

3. With probability qs we add m units of initial activ- 
ity on land that is previously undeveloped: the first 
end-point is selected without bias similarly to Eq. 
J[J but among undeveloped nodes. However, since 
nodes have no properties beside their degree in the 
non-spatial model, any undeveloped node can be 
added. The second end-point is selected preferen- 
tially according to Eq. 

We will refer to the growth classes as type-1, type-2 
and type-3 growth and to their relative rates as q±, q2 
and qs with qi + 92 + <Z3 = 1 throughout the paper. 

1. A continuum formulation 



A. Formulation of the non-spatial model 

A geographic area on which an urban economic system 
can grow is represented by an enumerated set of nodes, 
{1, 2, . . . , N} corresponding to non-overlapping land ar- 
eas. Trade of goods and service between activities in the 
nodes is represented by undirected edges. Since activities 
within the same node can trade with each other, an edge 
can connect a node with itself. The amount of activity 
of a site x t is defined as the degree of the corresponding- 
node. See figure ^ 

The network is initialized by connecting no nodes so 
that each has a degree of xq. The N — no undeveloped 
nodes have no trade interactions and thus no activity. 

At each time step we update the network as follows: 

1. With probability qi we add m edges between sites 
that already are developed: the first edge end-point 
is selected uniformly among developed nodes. The 
probability of a node i to be selected is 



Instead of assuming any particular set of activities and 
interactions we can use a continuum formulation where 
we consider each edge end-point to be an average of a 
large set of urban activities. It follows that an average 
activity must be assumed to interact equally much with 
all other average activities. Rather than counting explicit 
interactions to determine the activity level, we study the 
evolution of the expected node degrees. The time evolu- 
tion of activity on a developed site i follows the equation 



x l (t+l) = x l (t)+q 1 —r^+m(q 1 +2q 2 +q 3 )=-^ — -r, (3) 
n\ l^ l j x j\ t ) 

which is solved by the continuous-time method intro- 
duced by Barabasi et al [3. After sufficiently long time 
the degree distribution approaches the form 

P[ Xi = x ] ~ (x + A)-">, (4) 

with 



(1) 
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and 



According to Eq. J3J the node degree distribution will 
be power law distributed for the non-spatial model. 



B. Formulation of the spatial model 

The non-spatial model does not include any informa- 
tion about inhomogeneous relations in the network. An 
important departure from such a simple model that can 
be readily observed in real systems is of course the fact 
that different pairs of sites have different distances be- 
tween them. Because of transportation cost optimization 
this has the potential to affect the edge distribution of 
the network model. 

We will now extend the basic model to incorporate spa- 
tial interactions. This allows us to verify that the power 
law distribution of activities that is predicted by the non- 
spatial model is retained when space is introduced. Fur- 
thermore, it allows us to better map model output to em- 
pirical measurements. Scale-free spatial network models 
that have been studied recently are not directly appli- 
cable to urban economics since they require an a priori 
distribution of node locations HUGH 



1. Network evolution 

In the continuum non-spatial model we were able to 
treat network evolution as a purely local phenomenon 
since interactions in the network are homogenous, see 
Eq. © . In the spatial model we have to again explicitly 
separate the dynamics of the edge end-points. This is 
because the spatial context of the first end-point of a new 
edge modifies the probabilities with which other sites will 
become the second end-point (transportation costs). 

The same types of growth that are used in the non- 
spatial model are also used in the spatial model. But, 
as stated, the two end-points here need to be dealt with 
separately. The first site (primary increase) is selected 
either uniformly or preferentially depending on growth 
type. This will be discussed shortly. The second end- 
points of the trade relation edges (secondary increase) 
are always preferential. 



2. Secondary increase 

The spatial model tells us how a given primary increase 
Pj in activity at site j causes secondary increases Si in 
activity at all other sites i where 

s i=Pi v n _ . ( 7 ) 



This is analogous to Eq. J5J but in the spatial model, 
transportation costs will bias the choice of trade partners. 
Dij is a matrix representing the interaction strength, 
which is assumed to decay with the transportation costs. 
For the results in this paper we have used = (1 + 
cd(i,j))~ a , where d(i,j) is the Euclidean distance be- 
tween sites i and j. The non- negative parameters c and 
a controls the impact of spatial. Another function that 
can apply is = max(0, 1 — d(i,j)) if the transporta- 
tion characteristics of the activity is known, which can 
be the case in a model where activity types are mod- 
elled separately rather than as an average. Exponential 
= e~ d ( IJ ) can also be interesting to the extent that 
shielding is important, i.e. an activity tends to trade 
exclusively with the nearest supplier. 



3. Primary increase 

When trade takes place there is a mutual benefit in 
efficiency that is often used for further increasing the ac- 
tivity in the city. This feedback process makes it possible 
for the amount of activity in the urban system to increase 
considerably faster than the population. It is useful to 
think of Eq. J7J as a black box system to which a driv- 
ing force, primary increase, is applied. It should also 
be noted that the primary effects used for the model we 
present here are by no means neither exhaustive nor final: 
most earlier models of urban growth could be introduced 
as primary effects in our framework. 



4- Primary effects in type-1, type-2 and type-3 growth 

Network evolution in the spatial model is similar to 
the non-spatial model. Secondary increases are always 
preferential following Eq. 10 and primary increases now 
reflect the spatial distribution: 

Type-1 growth: the primary uniform increase is iden- 
tical to the non-spatial case following Eq. QJ. 

Type-2 growth: the primary preferential increase is 
identical to the non-spatial case following Eq. J2J). 

Type-3 growth is split into two related processes where 
one corresponds to growth in the perimeter of clusters 
and the other corresponds to growth in connection to 
infrastructure in the rural areas between clusters. Such 
infrastructure is not explicitly represented in our model 
and instead we use a parameter e to tune the amount 
of ambient infrastructure and thus the rate with which 
seemingly isolated clusters will appear. 

Type-3a growth: with a probability of 93(1 — e) we 
set the activity of a perimeter node to m. Perimeter 
nodes are nodes that are not developed but borders to 
at least one developed cell. The site of the new node is 
selected randomly and with uniform probability among 
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the perimeter sites on the grid 



(8) 



The total growth for a developed node can in this case 
be separated into a uniform and a preferential part as 



Xi(t + 1) = Xi(t) + Ci{t) + Vi(t)xi(t), 



(17) 



where nf is the probability with which node i is selected 



(P) 



1 if the node i is on 



to undergo type-3a growth, 5, 
( p) 

the perimeter, 6} = if the node is not on the perimeter 

(p) 

and n t is the number of perimeter nodes at time t. 

Type-3b growth: With a probability of q^e we set the 
activity of an external node to m. An external node is 
a node that is undeveloped and that has no developed 
neighbors. The site of the new node is selected randomly 
and with uniform probability among external sites on the 
grid 



1 



AE)> 



(9) 



where 7q is the probability with which node % is selected 
to undergo type-3b growth, 5\ = 1 if the node i is 

( E) ( E) 

external, 5} = if the node is not external and n t is 
the number of external nodes at time t. 



5. A continuum formulation of the spatial model 

A continuum formulation for the evolution of devel- 
oped nodes in the spatial model can, as in section III A II 
be constructed by studying the time evolution of ex- 
pected node degrees. 



Xi(t + 1) = Xi{t) + E\pi{t)] + E[ Si {t)}, 



(10) 



where E[si] is the expected secondary increase which can 
be calculated by a weighted summation over the expected 
primary increments, 



(11) 



For the special case of our simple model for the primary 
effects we have 



with 



A 2 )\ 



(P) (3a) . AE)(3b) 



6 i 'Pj 



P, 
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9i- 



,(D)' 



(2) mx 

Pi = 92 



,(3a) 



93(1 



, m 



AP) ' 



P 



(3b) _ 



93£" 



n 



(E)' 



(12) 

(13) 
(14) 
(15) 
(16) 



with 



and 



Ci = 91 " 



92 m 
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E[ Pj \ 



D: 



Efc D kj x k ' 



(18) 



(19) 



Noting that J2iVi x i 
be rewritten as 



m(qi + 2q2 + (ft), eauationll7lcan 



Xi (t + 1) = Xi(t) + (i(t) +m{ qi + 2q 2 + q 3 ) 



T]i(t)Xi(t) 

(20) 

which, in a comparison with equation J3J), reveals that 
the only difference between the the spatial and the non- 
spatial model is the site and time-dependent parameter 
rji. This is similar to the concept of node fitness, as pre- 
sented in 0, ^2 > which can affect the node degree dis- 
tribution. However, our simulation results indicate that 
rji falls within a sufficiently narrow interval for the power 
law to be essentially preserved (Fig. |5J). This is also 
supported by calculations of rji for both simulated and 
empirical data. 



C. The network model in an urban economics 
context 

1. The connection between node degree and land value 

An approximate linear relationship between node de- 
gree in the model and land value in the real system is 
crucial for the interpretation of our results. The motiva- 
tion follows from i) market pricing of goods and services 
and ii) the connection between trade benefits and land 
value. 

i) Market pricing of commodities provides an adap- 
tive measure that allows us to compare the activities 
that generate them. Hence, on average, an edge con- 
tributes identically to the value of both nodes to which 
it connects. This contribution is exactly our definition 
of activity, which implies that the degree of a node is 
proportional to its benefits due to trade. 

ii) This connection consists of two proportionalities. 
For a node i we have 



Vi OC n OC Xi, 



(21) 



where i>j is the value of the corresponding land area, r.; is 
the bid-rent 0Q, and Xi is the total trade benefits as 
outlined above. Capitalizing periodic rent income from 
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the site i gives land value Vi = where i is the inter- 
est rateQ. The second proportionality is a weak form 
of the leftover principle from urban economics, which 
states that, in a competitive land market, rent equals 
the amount of money leftover after all expenses (except 
rent) are paid. This amount of money equals the sum of 
all trade benefits at the site. For our results it is suffi- 
cient that, on average, a certain proportion of each new 
unit of trade benefit goes to the landowner. 

Together, i) and ii) suggest an approximately linear 
relationship between node degrees and land values. 

2. Types of growth 

Urban activity can increase in essentially two ways: ei- 
ther a new activity is related to, or it is unrelated to, an 
existing activity at the site. In the former case (pref- 
erential growth) this could be a new employee hired as 
a response to increased demand, in the latter case (uni- 
form) it could be the establishment of a new firm. Pref- 
erential growth corresponds to a per-unit activity rate. 
Uniform growth corresponds to establishment among lots 
on a competitive land market where, for the average land 
use, we can not expect any lot to be more profitable than 
any other. 

3. Reasons for treating perimeter growth separately from 
internal growth 

The jagged perimeter of urban areas exposes large ar- 
eas of undeveloped land to urban infrastructure, thus 
making it attractive for urban land use. Because perime- 
ter land is in ample supply and currently has a low rev- 
enue, even land uses with a very low trade gain can be 
competitive. Many low-activity land uses in the outskirts 
of the urban area can likely just barely out-bid agricul- 
ture and would not be viable in competition with other 
urban land uses. Among high- activity land uses some 
have very specific demands on land improvements and 
can therefore not benefit from buying existing buildings 
inside the urban area. This creates a special case for 
perimeter land. Note that just like for type-1 growth, 
competition prevents prediction of where the next growth 
event will take place among the perimeter nodes, and no 
preferential growth is possible since there is no previous 
activity that can expand. 

4- Node fitness and growth 

If we only regard trading activities the only difference 
between two sites with identical activity is the value of 
their spatial context. Therefore, in the spatial model, 
secondary growth is not homogenous, see Eq. Q. Node 
fitness (Eq. I17f) can be viewed as a local interest rate that 
predicts the growth rate of activity investments made in 



that site. The result of this is non-trivial growth predic- 
tions since the expected amount of new local development 
does not become a simple fraction of current develop- 
ment. 

Most notably, the model predicts the emergence of ur- 
ban sub-centers. This is realized by examining the ex- 
pected secondary effects in the spatial model (Eq. Illfl . 
Apart from being proportional to the amount of present 
activity E[si] is also subject to site competition and will 
increase for nodes that have high activity in relation to 
their own neighborhoods. For each possible primary ef- 
fect in j, the node development Xi is weighted with the 
fraction between trade intensity between site j and i and 
the sum of the all other trading options for the primary 
effect under consideration. Thus, nearby high-intensity 
nodes will not necessarily benefit a small neighbor. 



III. RESULTS 
A. Land values are power law distributed 

From empirical data, land values in Sweden are demon- 
strated to follow a power-law distribution for the higher 
range of land values (see figure [21 ■ The sharp transition 
that can be observed around 60kSEK suggests that two 
truly different mechanisms generate the prices below and 
above this point. This is in agreement with the observa- 
tion that the pricing mechanism we suggest would apply 
only to trading activity. 

The data we have used is based on the land value 
component of market value estimations of about 2.9 mil- 
lion units of real estate in Sweden. The data was origi- 
nally compiled by the Swedish National Land Survey and 
coded by Sweden Statistics to geographical coordinates. 
The data points we use are aggregated land values into 
IOOto x 100m squares. 

Our results are supported by a recent study by Kaizoji 
that shows scale-free behavior of land prices in Japan, 
with an exponent 7 ranging from 2.53 to 2.76|15|. 



B. Power law distributed prices are predicted by 
the non-spatial model 

We have developed a simple network model of the 
urban economy based on the Barabasi-Albert model 
by mapping fundamental assumptions from urban eco- 
nomics to the ontology of the network model (see section 
III A|) . In our derivation, we have determined the expo- 
nent of the node degree distribution, and thus, per our 
definition, also the predicted land values, to follow equa- 
tion JHJ. To the extent that our interpretation of the 
underlaying dynamics is correct this demonstrates why 
urban land values can be expected to follow a power law 
and how the exponent may depend on parameters. 
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tended to incorporate spatial constraints. 

The applicability of the network paradigm to urban 
growth suggests that scale invariance in urban systems 



V V 




FIG. 3: Simulated configuration from the spatial model with 
the parameter values used in previous plots. 



FIG. 2: Double logarithmic histogram of simulation results 
and empirical data. Simulation results are are denoted with 
squares in the figure and they are a mean of the results of 
three runs of the spatial model. The exponent of the model 
output has been tuned to match the empirical data. As indi- 
cated in Eq. \Q) this is done by setting the relative proportions 
of the growth types, in this case qi = 0.1, q2 = 0.6, e = 0.01, 
c = 0.2, a = I, m = WOkSEK and t =170000. The exponent 
is roughly 2.1 which is close to the value of 7 = 2.25 that 
is predicted by Eq. 0) for these parameter values in the non- 
spatial model. A slightly different value for the spatial model 
must be expected because of spatial bias in the growth dynam- 
ics. The sharp transition that occurs around a land price of 
OOkSEK (» 5kUSD) per 100m x 100m marks the difference in 
dynamics between trade based urban activities and rural activ- 
ities whose values are not described by the presented model. 



C. The spatial model retains power law statistics 



As discussed in section III B 51 the impact of spatial 
constraints closely resembles that of node fitness|llLll2j. 
This could potentially result in the distribution for the 
spatial network model becoming a sum of many power 
laws with different exponents. 

In figure [21 we plot results from simulations showing 
that node degrees in the spatial model follow a power law 
distribution. The model parameters have been tuned (see 
sec. Ill B 5|) to reproduce the distribution of the observed 
land prices. 



is caused by: i) growth and ii) preferentiality in how 
new trade connections are formed between areas of land. 
Growth in this context refers to the continual develop- 
ment of new land. Preferentiality is a consequence of 
point-to-point interactions between activities occupying 
the land areas. Note that many other observables, such 
as population and urban land use intensity, might be 
highly correlated with land value. 

The spatial model that we present can have more 
general applicability beyond urban economics. Other 
spatially growing networks are communication networks, 
transportation networks, electricity and utility networks. 
These can be expected to follow a similar type of growth 
since they are intimately connected to urban activity. 

The network architecture is generic and allows for ad- 
dition of any amount of detail. Also, being based on 
trade relations, the model produces output in units of 
currency. Because of this, such network models can pro- 
vide a bridge between a microscopic dynamics that can 
be found empirically and emergent economic properties. 

Further possible directions for research on urban eco- 
nomic networks include interpretation of other theoreti- 
cal network results in terms of urban dynamics, finding 
empirical parameters for scenario predictions and model 
validation. 
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IV. CONCLUSION 

We present empirical evidence that land prices follow a 
power law distribution for urban land uses, and present a 
generic model, based on the underlaying trade network, 
that reproduces this behavior. The model is a version 
of the Barabasi- Albert scale free network that is also ex- 
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